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Abstract 

Leptoproduction of nucleons into the backward hemisphere on nuclear targets is studied 
at relativistic subasymptotic energies and momenta. Spins are neglected. The relativistic 
internucleon potential is extracted from the appropriate photoproduction data. Different 
production mechanisms are shown to work together and interfere. Calculations show that 
whenever rescattering is possible it gives the bulk of the contribution except at very high Q 2 . 
The Weizsaecker- Williams approximation is found to generally reproduce only a small part of 
the total cross-section. Comparison with the data for A(e,e'p) reaction at E=2.4 GeV shows 
a reasonable agreement. 

1 Introduction 

Production of particles off nuclei in the cumulative domain x > 1 presents a great interest, 
since it is related to the nuclear structure at small distances, where one may expect formation 
of multiquark states and even droplets of a cold quark- gluon plasma. Staying within the most 
conservative picture of the nucleus consisting of nucleons interacting with some potentials, 
cumulative particle production allows to study the high-momentum asymptotics of these 
potentials and, which is of especial importance, multi-nucleon forces, otherwise hidden due 
to the low probability for several nucleons to be located at the same point. Experimental 
data demonstrate that the rate of cumulative production is very weakly dependent on the 
initial energy, starting from several GeV's. This allows to study cumulative production at 
comparatively low energy facilities, provided their luminosity is high enough to observe cross- 
sections dramatically falling as x grows beyond the threshold x = 1. 

There has never been shortage of theoretical models to describe cumulative production 
off nuclei. All proposed models can be roughly divided into three categories. The first (also 
in time) is to ascribe the production to the presence of states with a nuclear density much 
higher than on the average in the initial nucleus [1,2]. Interpretation of these states may 
be different: just several nucleons at the same point ("many-body correlations" [3] ) or multi- 
quark states [4]). The second mechanism assumes that such high-density states are rather 
formed in the course of the collision ("compressed tube" models [5]). Finally one can assume 
that the high nuclear densities are irrelevant and the cumulative particles are produced as a 
result of simple rescattering [6] . 

Microscopic studies, both in the nucleon and quark pictures, have shown that in reality all 
three mechanisms work together in the cumulative production. Moreover their contributions 
interfere, so that the total rate is not just the sum of them [ 7 ]. It is of most importance 
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that the contribution from the rescattering depends differently on the target properties that 
the other two. This prevents simple estimates of the cross-sections on heavy nuclei based 
on those with only few nucleons, as suggested in [ 8 ]. Only at sufficiently high momentum 
transferred to the nucleus, of all contributions, only the first one remains, corresponding to 
the few-nucleon correlations or multi-quark states inside the nucleus. From the theoretical 
point of view, it is the most interesting part, since it reveals the inner structure of the nucleus 
at short distances. However to see it clearly against the background of other components, 
one has to select events with high enough Q 2 , which drastically diminishes already small 
cross-sections. Also the exact values of Q 2 at which one can neglect the contributions from 
other two mechanisms are not known apriori. 

For this reason it is important to carry out calculations of cumulative production at 
finite Q 2 which take into account all three above mentioned mechanisms. In this way one 
expects to obtain predictions for the inclusive cross-sections, integrated over all Q 2 , and thus 
directly related to the bulk of experimental data. Moreover, studying the relative weight of 
the contributions at different Q 2 one may determine the minimum Q 2 value, starting from 
which rescattering and compressed tube contributions can be safely neglected. 

In this paper we do these calculations in the framework of the nucleon degrees of freedom. 
Estimates have shown that this description gives reasonable results at x not too close to 
the cumulative threshold x = n for n nucleons which interact at close distance from each 
other, provided one uses the correct relativistic kinematics [ 3 ]. We limit ourselves with the 
cumulative region I < x < 2, easiest from both experimental and theoretical points of view. 
In this region it is commonly assumed that one has to study hard interactions of only two 
nucleons inside the nucleus. Since inclusion of relativistic spins presents serious difficulties for 
nuclei with A > 2, 3, we simplify our approach by taking nucleons scalar. This approximations 
may certainly change results by a factor of order 2-^3. However in the cumulative region as 
x grows from 1 to 2 the cross-sections fall by several orders of magnitude. On this scale a 
factor of order unity plays a secondary role. In any case there is no reason to expect that 
the results obtained with scalar nucleons will have a different x, p\ and A dependence as 
compared to the full spinor case, and it is this dependence which is of primary theoretical 
interest. 

The paper is organized as follows. In Section 2 we present our basic formulas for the 
structure function and solve the most difficult task of extracting the rescattering contribution. 
In Section 3 we discuss how to pass from the structure function to the production rates. Our 
final results are presented in Section 4. Section 5 contains some discussion and conclusions. 
Appendix is devoted to certain non-trivial details of numerical calculations. 



2 Structure functions in the interval 1 < x < 2 

2.1 Cross-sections and kinematics 

We shall be interested in the inclusive leptoproduction of nucleons on nuclear targets in the 
cumulative kinematical domain. However it will be convenient to start with a simpler reaction 

e + A->e + X. (I) 

Its cross-section can be expressed via nuclear structure functions in the standard manner 

VMott^rj (^=F 2 {x, Q 2 ) + 2t&n 2 ?-F 1 (x,Q 2 ) ) , (2) 



dE'dVL' AttM \qP v ' 2 

where 



*Mott = ^(E>) 2 cos 2 -. (3) 
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In our normalization the structure functions are related to the imaginary part of the forward 
amplitude for the elastic 7* — A scattering as 



WW = {-9r + q ~f) F ^Q 2 ) + ^ - *rp) {P» - F 2 (x,Q 2 ), (4) 

where q and P = Ap are the momenta of the photon and nucleus and, as usual, 

QW, .-g. (5) 

Note that we define x respective to 7* + N scattering, so that 

< x < A (6) 

The region x > 1 is cumulative. 

The standard way to find the structure functions from W^ v is to choose a coordinate 
system in which q + = q y = pj_ = ("theoretical"). Marking the components of vectors and 
tensors in this system with bars, one finds 



0^ 

2Axp'i 



F 1 (x,Q 2 ) = W yy , F 2 (x,Q 2 ) = w ^ IT W ++ . (7) 



We take the nucleon at rest: p = 0. To study leptoproduction of cumulative nucleons we 
shall have to use another system ("lab") in which the incoming lepton with momentum I is 
moving along the z-axis in the opposite direction, so that neglecting its mass = 0. These 
two systems are related by a spatial rotation, which can be split into two ones. Starting from 
the lab. system, the first rotation is to put the y-component of q (or 1') to zero. This is just 
a rotation in the xy plane by the azimuthal angle <p: 

q' x = q x cos <f> + q y sin cf), q' y = -q x sincp + q y cos <\> = 0, q' z = q z . (8) 

The second rotation in x z plane is to put q + = 0, that is q z = — q^: 

q" = q' z cos <pi + q' x sin fa = —q , q x = -q' z sin fa + q' x cos fa . (9) 

The first relation determines the value of angle fa through the given value of the transferred 
energy qo = E — E'\ 

cos fa = —(-qoq z ± Qq' x )- (10) 

The two signs correspond to the two possibilities of directing q x parallel or antiparellel to the 
x-axis. 

To find the inclusive cross-section for the leptoproduction of cumulative nucleons we shall 
integrate the cross-section (2) over E' and O , present both structure functions in (2) as 
integrals over intermediate particles and then lift the integration over the momenta of the 
observed nucleon. The only point to have in mind is that the structure function and the 
inclusive cross-section are to be calculated in different reference systems, theoretical and lab 
respectively. Since both Q 2 and qp will vary in the resulting expression for the production 
rate, it will be more convenient to characterize the cumulative region directly in terms of the 
momentum p\ of the observed nucleon in the lab system. We shall discuss the corresponding 
kinematics in Section 3. 

The method we use to find the cross-section for the inclusive leptoproduction of nucleons 
from the inclusive leptoproduction of leptons is discussed and justified in Appendix 1. There 
one can also find the expression for the double inclusive cross-section for the process in which 
both the lepton and cumulative nucleon are observed. This cross-section can be found from 
(2) only averaged over the azimuthal directions of the observed nucleon. 
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2.2 Hard interactions between two nucleons 

In the interval n — 1 < x < n the incoming virtual photon has to interact with at least n 
nucleons. Since the cross-section results rapidly diminishing with the number of interacting 
nucleons, our first approximation will be to take into account only the interaction with n 
nucleons in this interval. Then at 1 < x < 2 the incoming photon is to interact with two 
and only two nucleons (approximation of "pair correlations" for the short-range part). The 
two interacting nucleons will generally have large momenta k and q — k, wheras the rest 
intermediate particles will be slow. Their momenta will be of the order of the standard 
nuclear ones, to be neglected wherever possible. Separating the large momentum part from 
the low momentum one we can present the total contribution to the ampliude for the reaction 
7* + A — ► N(pi) + N{p 2 ) + X as a sum of four diagrams, shown in Fig. 1, a — d, where the thin 
lines represent slow nucleons, thick ones represent fast or highly virtual nucleons and thick 
dashed lines represent the interaction between the nucleons with high momentum transfer. 

An important point is that the high-momentum part responsible for the production of 
the cumulative nucleons involves both short- and long-range distances in the nucleus, the 
latter corresponding to the rescattering contribution. Also to note is that the short-range 
contribution comes both from diagrams a and b and from diagrams c and d. The first pair of 
diagrams corresponds to the presence of the correlated pair inside the initial nucleus before 
the collision. Ideologically this is close to the "flucton" model, which was used to explain the 
cumulative phenomena immediately after their observation. The second pair of diagrams, c 
and d, correspond however to correlations created in the course of the collision. This is in 
line with mechanisms proposed in so-called "hot' models, like the "compressed tube model". 
As we observe, inspection of Feynman diagrams shows that in fact all different mechanisms 
proposed so far to explain the cumulative phenomena coexist simultaneously. Moreover, as 
we shall discover presently, their contributions to the cross-section are not added incoherently, 
but interfere, which leads to important cancellations. Our task will be to separate the long- 
range rescattering contribution, and to take into account the mentioned interference, which 
we shall do following [ 6 ] . 

The imaginary part of the forward scattering amplitude corresponding to Fig. 1 is given 



(2vr) 4 5 4 (fe 1 + k 2 - k[ - k' 2 )(2ir) 4 5( Pl +p 2 -q-k 1 - k^H^HlyT^MKA), ( n ) 



where the notations are clear from Fig. 1. In particular f(Q 2 ) is the effective nucleon form- 
factor; blob T which includes the propagators of the nucleons 1 and 2 describes the nuclear 
recoil particles. 

Factors H and H' depending on kj and kj, j = 1,2 respectively, correspond to the high- 
momentum part of the amplitude. Their denominators in the virtual nucleon propagators in 
Fig. 1 a and b are respectively 



Obviously they cannot vanish, since a stable nucleus cannot decay into a real nucleon and 
some other particles. The denominator in Fig. 1 c and d is 



by 




m 2 — (2p — p\) 2 , and m 2 — (p\ — q) 2 . 



m 



2 -(q + h) 2 = [S(ki + q) - E(k 2 ) + E( Pl ) + E(q - Pl + kx + k 2 ) 
[E(ki + q) + E{k 2 ) - E{ Pl ) - E(q_ - P i + kx + k 2 )" . 



(12) 
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The first factor is the difference between the initial and final energies in the decay 

7V(k 2 ) - 7V(-kx - q) + AT(pi) + iV(q - Pl + k x + k 2 ) 

and cannot vanish due to stability of the nucleon. For this reason we can drop small momenta 
ki t 2 in it. The second factor is the difference between the initial and final energies in the 
rescattenng 

" iV(ki + q) + iV(k 2 ) - 7V( Pl ) + AT(q - Pl + k x + k 2 ) 

and evidently may vanish. Taking into account the smallness of momenta k\ and k 2 we can 
rewrite it as 

cv \ , EV a it 1 / \ k iq (ki+k 2 ,q-pi) . . 

£(q) + m - £( Pl ) - E(q - Pl ) + _ - ^— . (13) 

So finally the denominator (12) takes the form 

m 2 - (q + ki) 2 = a(u + a - iO), (14) 

where 

a = £(q)+£(pi)+£(q-pi)-m, (15) 
u = £(q)+m-£(pi)-£(q-pi) (16) 

and 

= kiq _ (ki + k 2 ,q-pi) 
° £?(q) £7(q-pi) ' 1 ' 

Thus approximating k\ = k 2 = p everywhere except in the denominator (14) we find the 
high-momentum factor as 

^ m 2 -(2p-pi) 2 m 2 -(q-p 1 ) 2 a(u + a-iO) v 7 

We are only interested in + and y components of H. Factors K + and K y are identical for H 
and H' and in our system they are 

=2(2p + -pi+), =2pi+, K| c) = 2p+ (19) 

and 

#W = = -2p ly , = . (20) 

To move further we first integrate over the zero components of fcj and fcj. The (5 functions 
which correspond to the conservation of energy have the form 

S(k w + k 20 + P^- P )S(k 10 + k 20 + P^- P )S(p 10 + P20 + Po~qo- P ). (21) 

Here Po and Pq are the energies of the target nucleus and its debris respectively. Using the 
first two (5-function we integrate over the zero components of the two active nucleons 1 and 2. 
In coordinate space this integration puts both active nucleons at the same time. If we take 
into account only poles of their propagators and free nucleons for the debri, it produces the 
standard denominators of the wave function of the target nucleus. 
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Now we make the decisive approximation and neglect small energies of the nuclear debri 
as compared to fast particles in the last (5-function in (21) taking it as 

<5(Pio + P20 -qo~ 2m). (22) 

This allows to integrate blob T over the energy Pq of the debri, which puts initial and final 
active nucleons at the same time point. As a result T transforms into the nuclear /j-matrix 
for a pair of active nucleons: 



/ 



dk20dk> 20 dP{> T{k ^ k ^ k ^ = _ 1)p(ki)k2 | k / )k / 2) _ (23) 



2vr 2vr 2vr 

It is convenient to pass to the p matrix in the coordinate space presenting: 
(2nf5 3 (ki + k 2 - k[ - ^)p(ki, k 2 |ki, k' 2 ) 

2 



/ 



J] d 3 r j (Pr' j p(T 1 , r 2 |ri, r^yi^+k^k'^y^) ^ 

3=1 

and write the imaginary part of the forward scattering amplitude in the form 



W, v = l -A{A-l)f(Q")JdrjY{ 



lr dPkjdPktjcPrjcPr'j 



/9 (r 1 ,r 2 |r / 1 ,r 2 )^ kiri+k2r2 - k i r i- k ^)F M (i70*, (25) 

where 

dT = (o^% (2vr) 4 ^i + P 2-q-2p) (26) 

(27r)- i 2pi (27r)- i 2p 2 o 

is the standard phase volume for the virtual photo-production on a pair of nucleons at rest. 
2.3 Short and long range parts 

We start with the component W yy . As follows from (20) it does not contain any terms 
singular at kj = kj = 0, j = 1, 2 and so can be taken at this point. Then integration over kj 
and kj, j = 1,2 gives four 3-dimensional S- functions which will put coordinates Tj = r'j = 0, 
J 1-2: 

W yy = l -A{A - l)f(Q 2 ) Po J drR 2 y , p = p(0, 0|0, 0). (27) 

Obviously this corresponds to a pure short range mechanism, both direct and spectator. 
Passing to W++ we split H + into parts singular and non-singular at kj = 0, j = 1, 2 

H + = H + + , (28) 

where H + is a non-singular part and according to (18) 

R = 9n ^(Pl -p)+v(pi-q-p) , , 

P+ E(q)+E(p 1 ) + E( (l -p 1 )-m- { > 
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Doing the same for (H' + )* we find 

H + (H'Y = Hl + 2BH + Re — !— + B 2 - , (30) 

+ v +; + + u-iO u + a-iOu + a' + i0' v ' 

where we have put kj = kj = 0, j = 1, 2, wherever possible. 

The first two terms in (30) do not depend on kj, kj, j = 1,2 and so integration over these 
leads to a short-range contribution similar to (27) 

W ++ = ±A(A-l)f 2 (Q 2 )p J dr^H 2 + + 2BH + Re-^y (31) 

To treat the singular contribution we first symmetrize it in primed and non-primed mo- 
menta 

x s 1 ' . ~l( 1 . 1 , + 1 



u + a — iOu + a 1 + iO 2 \u + a — i0u + a' + i0 u + a + iQu + a' — iO 
Next we use 

1 1 

+ 2tt5(u + a) 



to find 



X 



u + a — iO u + a + iO 

= - ( \ h 2iri5(u + a) \ 

2\u + a + i0u + a f + t0 u + a' + iO 

1 1 1 

H - 2mo{u + a 



w + a — iO u + a' — iO u + a' — iO 

= Re- ^ + 2vr 2 (5(u + a)5(u + a / )- (32) 

(u + a-i0)(u + a'-i0) v ; v ; v > 

The first term is non-singular at kj = kj = 0, j = 1, 2 and can be rewritten as Re l/(u — 
iO) 2 . The second reduces to 2ir 2 5(u)5(a — a'). Taking into account that k\ + k 2 = k[ + k' 2 we 
find that a — a' = (q, ki — k' 1 )/E(q), so that finally 

X = < 1 -^2 + 2vr 2 ^(q)J(n)5(q, k x - k' x ). (33) 
(u — 

The first term in (33) does not depend on kj = kj, j = 1,2 and its integration is similar 
to the previous terms. It also corresponds to a short-range contribution and summing with 
(31) forms the complete short-range part of W++: 



WH = ^A{A - 1)/ 2 (Q> / dr [h% + 2BH + Re — + £ 2 Re j-± 



(34) 



The second term in (33) describes the long-range mechanism corresponding to real rescat- 
tering. It contains non-trivial dependence on small nuclear momenta. To do the integrations 
over kj = kj, j = 1,2 we direct q along the z-axis. Then integrations over k±±, k' 1± , k2, k' 2 
and the corresponding coordinates are done as before and the only non-trivial integrations 
left are over k\ z and k' lz . We find the long-range part of W ++ as 

wH = \a{A-1) I dr27T 2 ^B 2 5(u) f ^lldz 1 dz' l e tk ^- z ^ p(0 ± , z V , 0|0j_, z[; 0) 
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2,£(q) 



where 



= 2vr Z-j^p y dTB 2 5(u), (35) 

/ = - 1) I d^(0 ±) z; 0|0j_, z; 0) (36) 

has a meaning of the average dimension of the nucleus, travelled by the nucleon in rescattering. 
Note that due to energy conservation u and a can be rewritten as 

u = E(q) -q -m, a = E(q) + q + m (37) 

and do not depend on integration variables in (35). So one can take 5(u) out of the integral 
(35). This leads to our final formula for the long-range contribution 



WH = 2K 2 lf 2 {Q 2 )5{E{q) - g - f drB 2 . (38) 

Note that it has a 5-like behaviour in x. Indeed 



6{E(q) - qo - m )= 2(q + m)5(E 2 (q) - (q + mf ) 

Q 2 + 2m 2 



= 2{q + m)5(Q z - 2mq ) = 5(x - 1) 



2mQ 2 



So in fact the rescattering does not give any contribution to the structure function in the 
cumulative region. However it does contribute to cumulative particle production. 

Note finally that our nuclear p- matrix is defined in the relativistic normalization. Its 
relation to the non-relativistic one is given by 

P = ^PNR. 



3 Cumulative nucleon production 

We shall be interested in the observation of a nucleon of momentum p\ produced in the 
reaction 

l + A^N( Pl )+X, (39) 

where I is the incoming lepton of momentum I and A is the target nucleus at rest. We 
shall direct 1 antiparallel to the z-axis to have the system close to the theoretical system in 
which q + = 0. With such a choice non-cumulative nucleons will be produced in the backward 
hemisphere, whereas cumulative ones will be observed also in the forward hemisphere. The 
exact non-cumulative kinematical region in the lab. system can be easily established in light- 
cone variables, using the fact that in this system l + = 0. Let p\ + = x\p + . Obviously in the 
non-cumulative region < x\ < 1. Conservation of the "-" component of the momentum 
(light-cone energy) requires 

1-+P- = Pi- +P-, (40) 

where p' is the momentum of the undetected particles. The minimal mass of the latter is 
zero. So in terms of initial energy E we get a relation 

'(2E + m)>^±^ + — ^ , (41) 
V2 2x lP+ 2(1 -xi)p+' 
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which leads to the condition 

pj± <xi(l- x 1 )m(2E + m) - (1 - xi)m 2 (42) 
Together with the relation between x±,p± and p\ z : 



Eq. (42) determines the non-cumulative region in the p\ z ,Pi± plane. The cumulative region 
for the pair correlation is determined by the corresponding energy conservation relation 

Z_ + 2p_ =p_ +p'_, (44) 

where now the minimal mass of undetected particles equals m. In terms of x\ and pi± this 
gives a relation 

p 1± = xi(2 — xi)m(E + m) — m 2 , (45) 

which together with (43) determines the cumulative region for the pair correlation. For the 
incident energy E = 6 GeV in the backward hemisphere this region is limited by the solid 
line in Fig. 2 The dotted line limits the region corresponding to the rescattering discussed in 
the previous section. 

As mentioned shall find the rate of cumulative nucleon production by using our results for 
the structure functions. Putting them into the integrated cross-section (2) we then lift the 
integration over the momentum p\ of the intermediate real nucleon. However we have to take 
into account that our expressions for the structure functions have been obtained in a special 
coordinate system q + = q y = different from the lab. system for the reaction (1). This 
implies that we have to reexpress all momenta in the system q + = q y = via the momenta 
in the lab. system. 

Turning to our expressions for the structure functions, we see that in most places the 
vectors appear in the form of scalar products or invariant integration measures, which do 
not change under rotations. The only quantities which change their form under rotations are 
vectors K, in which we shall have to express the +- and y-components of vector p\ in the 
system q + = q y = by its components in the lab. system: 

Ply -> Ply = -pix sin 4>+pi y cos (f>, p\ z -> pi z = p\ z cos 0i + (pi x cos (f) -\- piy sin <p) sin <f>i (46) 

and keeping the same po- 

Thus we get for the inclusive cross-section to observe a nucleon of momentum p\ 

<H^1^ = l^na MM ^(^F 2iq , Pl)+ 2 t ^F M , Pl) y ,47) 

Here the structure functions iq,2(q, Pi) are obtained from the standard ones, calculated in 
the preceding section, by lifting the integration over pi, that is discarding d 3 pi/((27r) 3 2pio in 
the integration volume dV in (34) and (38) , and also changing p\ y and p\ z to p\ y and p\ z as 
indicated in (46). Note that angles 4> and 4>i in the latter formulas depend on the remaining 
integration variables, related to the changing vector q = 1 — 1'. 
Using (7) we rewrite (2) in terms of W's to obtain 



(27r) 3 2pio<icr 1 



j dE'dna M ott (^W ++ {q, Pl )) + 2tan 2 ^W yy (q, Pl ) j . (48) 



d 3 pi Air Am J \P+ 
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In W's we have to lift the integration over pi. We rewrite the phase volume (26) in the form 
and integrating over P2 find 



j3 

dr = 3 Pl 2tt5((2 P + q - Pl ) 2 - m 2 ), (49) 

so that lifting the integration over pi substitutes 

dr -^2ir5((2p + q-p 1 ) 2 -m 2 ). (50) 

In variables E', 9, cf> we have 

qo = E — E', q z = —E + E'cos9, q x = — E'sm9cos(f), q y = — E'sm9sm4>, 

Q 2 = 2EE'{1 - cos 9). (51) 
The fixed momentum pi can always be placed in the xz plane, so that p\ y = 0. The argument 
of the (5-function in (49) turns out to be 

(2p + q - Pl ) 2 - m 2 = Am 2 - Q 2 + 4pq - 4p Pl - 2q Pl == 2(6 - E'\(9, (f))), (52) 

where 

A = E(l — cos 9) + 2m — E\ — p\ z cos 9 + p\ x sin 9 cos <p (53) 

and 

b = 2m{m + E - E{) - E(E 1 +pu). (54) 



Here E\ = y m 2 + p\ z + p\ x is the energy of the observed nucleon. Integrating this function 
we fix E'(9,(f>) and obtain an extra factor 1/(2A). Thus we get 

(27r) 3 2p 10 da 1 f d cos 6 'deb ( ' m 2 - , XN ,9- , A 

where it is understood that we have to drop dr in W's, make the substitutions (46) and 
express E 1 ' in terms of 9 and using (52): 

E' = b/X(9,<p). (56) 

Both 6 and \(9,<p) result positive and do not vanish in the integration region. However 
&Mott contains denominator Q A typical for virtual-photon-mediated processes and leading to 
a logarithmic divergence for massless leptons mi = 0. To cure this divergence we take in (55) 
the expression for Q 2 for non-zero mf. 

Q 2 = 2EE' - 2^(E 2 -mf){E' 2 - mf) cos 9 - 2m 2 . (57) 

Note that the standard Weiszaecker- Williams approximation takes into account only the 
leading term in log(l/m/) and thus expresses the cross-section via the real photoproduction. 
However calculations show that the Weiszaecker- Williams contribution is generally quite a 
small part of the total cross-section due to the fact that the amplitude vanishes at Q 2 = 0. 
This implies that leptoproduction cross-sections cannot be directly reduced to photoproduc- 
tion ones. 

Lifting one or both integrations in (55) and introducing appropriate phase volume factors 
one obtains cross-sections with Q 2 or 4- vector q fixed respectively, averaged over the azimuthal 
directions of the observed nucleon. If one is interested also in this azimuthal dependence, one 
should use the formula for the double inclusive cross-section derived in Appendix 1. 
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4 Photoproduction of cumulative protons and the choice of 
the potentials 

Expressions for H contain the relativistic potential v(k) which desribes the pair interaction 
between the nucleons. A variety of relativistic potentials have been proposed which correctly 
describe the nucleon-nucleon scattering data at moderate energies. However, on the one hand, 
we actually need the high-energy asymptotics of the potential, which is fixed by these data 
not well enough. On the other hand, the proposed potentials refer to relativistic nucleons with 
spins and it is not clear how they can be used in our spinless picture. For this reason we shall 
recur to a different approach. We shall use the fact that for the process of the production of 
cumulative protons by real photons, the cross-sections are directly related to the internucleon 
potential. Therefore staying in the spinless picture and using the experimental data on the 
photoproduction on deterium one can extract an effective internucleon potential as a function 
of the momentum transfer in a straightforward manner. However the energy dependence of 
this potential cannot be determined from the data, since the extracted potential refers to 
the energies close to the threshold. We studied two possibilities: an energy independent 
potential, similar to the exchange of scalar mesons, and a potential which grows with energy 
like the exchange of vectors mesons. Note that from a purely theoretical point of view the 
first alternative seems preferable, since it leads to vanishing of finite state interactions at high 
Q 2 in correspondence with colour transparency following from the QCD. Also, as we shall see, 
our numerical results with the energy independent potential agree quite satisfactorily with 
the existing data on cumulative nucleon production, whereas with the vector-like potential 
they overshoot the data by factor 2.5 (see Fig. 11 below). For these reasons, we discuss below 
only the extraction from the photoproduction of an energy independent potential. 
For unpolarized photons the total cross-section is 



a tot 



pol 



where J = AApq is the invariant flux, e M is the photon polarization vector and W^ v are the 
same as in Section 2. 

We direct the momentum of the coming photon along the z axis. The integrand for 
W^v involves a product H^H'* so that summing over polarization we get (eH)(eH)* = 
H(H')* ± . From (18) we find that diagrams c and d do not contribute to the cross-section in 
our approximation. Assuming that at Q 2 = only the proton interacts with the incoming 
photon, we are left with diagram b only (the direct mechanism) for which 

H± = 2 Pll _ . 59 

2qpx 

This obviously implies that there is no rescattering and we can safely drop small nuclear 
momenta in H. Repeating our derivation in Section 2 we find the total cross-section as 

a tot = h A _ l)po Jl f drH\. (60) 



2 V /ru 8pg 

Passing to the inclusive cross-section for the production of a proton we make the substi- 
tution (50) to find 

= \(A - D»^( (2P + q - Pl f- „?)Hl. ,61) 
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We shall use the lab system where the nucleus is at rest. Then qp = mE where E is the 
photon energy. The argument of the 5- function in (61) is t + AmE — AmE\ + 3m 2 , where 

t = (q - Pl f = m 2 - 2E(E 1 -p 1 cos6). (62) 

Ei is the proton energy and 6 the angle between the proton and photon directions. 

Transforming the phase volume element to variables t, Ei and using the 5 function to 
integrate over E\ we find 

Note that in H± the argument of the potential is 

(pi - q - V? = t + m 2 + 2m{E - Ei) = (l/2)(/ - m 2 ) = (l/2)*i 
and the denominator is 2qpi = —t±, so that we finally obtain 

For the deuteron we have to substitute 

±A(A-l) Po = (l/m)ij 2 (0), (65) 

where ip(r) is the deuteron wave function in the coordinate space. 

As we observe, the photoproduction cross-section reveals a great simplicity in our ap- 
proach: due to kinematical reasons the rescattering is absent and the cross-section becomes 
directly related to the high-momentum asymptotics of the effective internucleon potential. 
Therefore it is most suitable for studying this asymptotics. Of course one has to remember 
that our formulas may pretend to be valid only for the values of x sufficiently higher than 
unity. This is achieved in the emission of protons in the backward hemisphere relative to 
the direction of the incoming photon. Present data mostly refer to forward directions. Only 
the data [ 9 ] on the deuteron at the incident photon energy 4.0 GeV and 6\ = 90° involve 
sufficiently high values of x. These data show that with more or less standard nuclear poten- 
tials borrowed from low energy physics one cannot desribe the experiment above E = 1 GeV. 
The experimental data for s u da/dt show a clear plato above 1.4 GeV in accordance with 
the QCD scaling law, whereas the predictions from the meson-exchange potentals steadily 
rise [ 9 ]. This is not at all unexpected, since one cannot hope that the relatively low energy 
parametrizations for the potential will work sufficiently well at high momentum transfers 
involved. 

As mentioned, we are going to use the relation between the photoproduction cross-section 
and potential in the opposite direction. We determine the asymtotics of the effective nuclear 
potential from Eq. (64) putting the experimental data from [ 9 ] into its left-hand side. The 
latter can be satisfactorily fit by the expression 

s 11 -^ (kbn * GeV 20 ) = 0.4 + 0.147741e _E2 + 1.47116e _2E2 (66) 

valid for E = 0.2 4- 4.0 GeV. We take the value for the deuteron wave function at the origin 

as 

V>(0) = 3.131 10 _2 GcV 3/2 . 

The momentum space potential extracted in this manner is shown in Fig. 3. Compared to 
standard internucleon potentials it falls more rapidly with the transferred momentum. This is 
in agreement with the reduced nuclear amplitude approach of [ 10 ], in which hard exchanges 
between nucleons are to be accompanied by damping form-factors. 
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5 Numerical results 

Using the effective potential extracted from the photoproduction data we calculated the 
leptoproduction cross-sections (55) for the initial lepton energy 6 GeV in the backward hemi- 
sphere on various nuclear targets. The A dependence of our cross-sections is concentrated in 
factors po and I for short- and long-range contributions correspondingly. These factors were 
calculated in the standard manner, neglecting correlations in the nuclear p-matrix and using 
the Woods-Saxon nuclear density. One finds 



Calculated dimensionless /o and Co are shown in Fig. 4 for different nuclei. 

Numerical integrations in (55) encounter certain difficulties related to the singularities of 
the integrand at Q 2 = and the values of angle 6 limiting the rescattering region. Some 
details on this point are given in the Appendix 2. 

Our results for the inclusive cross-section are presented in Figs. 5 and 6 for Cu and Pb 
targets for various angles in the backward hemisphere. 

To see the relative importance of different production mechanisms, in Figs. 7-10 we show 
differential cross-sections in both pi and Q 2 at Q 2 = 1 and 4 (GeV/c) 2 and 90° and 180° on 
the Pb target. Solid curves correspond to the total contribution. Dashed and dash-dotted 
curves show contributions from the spectator mechanism and rescattering correspondingly. 

As a general result we find that for the studied values of Q 2 in the kinematical region 
where rescattering is possible it either completely dominates (at Q 2 ~ 1 (GeV/c) 2 or lower) 
or gives a contribution of the same order as the spectator mechanism. The contribution of 
the direct mechanism has been always found quite small. At comparatively low values of Q 2 
the short range contribution from the final interaction ("compressed tube" mechanism) is 
also of the same order as the spectator mechamism and negative (see Fig. 8). However with 
the growth of Q 2 its relative weight diminishes, as expected (cf. Fig. 10). We have also found 
that the relative contribution from the Weizsaecker- Williams term (leading in log(l/m;)) is 
generally quite small, except at highest values of p\ for a given angle where it becomes of the 
same order or even several times larger than the rest short-range contribution. 

To have some contact with existing experimental data on A(e,e'p) reaction we also calcu- 
lated the double inclusive cross section da / {dE' dVLdVLi) on the 16 target, where fl and Qi 
are angular variables for the final lepton and nucleon respectively. We have taken E = 2.4 
GeV, Q 2 = 0.8 (GeV/c) 2 and qo = 0.439 GeV corresponding to the experimental kinematics 
in [ 11 ]. Our results for the backward hemisphere are shown in Fig. 11 together with experi- 
mental points for the protons knocked out of lpi/2 shell (lower points) and IP3/2 shell (upper 

points). The recoil nucleon momenta result small in this kinematics, so that our approxima- 
tion of neglecting the binding does not seem to be well justified. However the agreement is 
unexpectedly good: our calculated cross-sections well correspond to the experimental ones 
averaged over the two nuclear levels. 

6 Conclusions 

We have studied leptoproduction of cumulative nucleons on nuclear targets in a realistic 
subasymptotic kinematics, where all proposed mechanisms work together and interfere. Our 
results show that in the kinematical region where rescattering is possible, its contribution 
cannot be neglected unless Q 2 is fixed and large. For the inclusive cross-section intergrated 



(A - l)po = ^A(A -I) J d\p 2 NR {r) = m 2 f ( 



(67) 




(68) 
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over Q 2 rescattering dominates the cross-section in this region. Of the rest mechanisms the 
direct one has been found to be completely unimportant at Q 2 > 1 (GeV/c) 2 . However 
at smaller Q 2 its contribution interferes with the rest mechanisms to make the amplitude 
vanish at Q 2 = 0. The short range part of the final state interaction (the compressed tube 
mechanism) gives a negative contribution, of the same order that the spectator contribution, 
except at high values of Q 2 when it diminishes. 

So, on the one hand, our results confirm the standard expectations that at very high Q 2 
only the spectator mechanism survives, which allows to simply relate the cumulative produc- 
tion on different nuclear targets [ 8 ]. On the other hand, they show that for realistic not 
too large Q 2 and certainly for the inclusive cross-sections integrated over all Q 2 all mecha- 
nisms give contributions of comparable order and interfere. In particular, in the kinematical 
region accessible for rescattering the latter gives a sizable (or even dominating) part of the 
contribution, which prevents using simple scaling arguments to fix the A-dependence. 

Another important result is that the Weizsaecker- Williams approximation generally con- 
stitutes only a small part of the total cross-section. Thus leptoproduction cross-sections 
cannot be simply reduced to photoproduction ones. 

Our calculations have been based on a simplified picture, in which all spins have been ne- 
glected. Accordingly the effective internucleon relativistic potential have been taken from the 
experimental photoproduction cross-sections, calculated within the same spinless approach. 
Comparison with the existing data on A(e,e'p) reaction reveals a reasonable agreement, which 
supports validity of our treatment. 
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8 Appendix 1. The double inclusive cross-section 

In principle lifting integrations in the inclusive cross-sections does not give the exclusive ones: 
one can add any function to the thus obtained cross-section which integrates out to zero in 
the inclusive cross-section. So to have the double cross-section for the process in whichboth 
the lepton and nucleon are observed we have to study the amplitudes as they are, not via 
the inclusive form with the structure functions F\^- In our case this is possible (since the 
amplitudes are known). 

The starting point is the basic formula for the cross-section 



do = -J^fL.irWrdr (69) 



in which API is the flux, L and W are the leptonic and hadronic tensors and dr is the phase 
volume. The leptonic tensor is 

V = \ E^VX^V) = 2(^„ + IX - (n')g^). (70) 

pol 

The hadronic tensor can be represented by a symbolic product of two currents J which are in 
fact matrix elements of electronagnetic current between the initial and final hadronic states: 

wv = <v;. (7i) 
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Thus the cross-sections has the form 

1 e 4 



da = 



2(Pl) d 



-dr 



(i'j)(uy + (u)(i'jy - {u'){jj*) 



(72) 



We choose the system l + = l± = 0. Then both the numerator and denominator become 
proportional to /_, which cancels out. We get 



da 



1 e 4 



d,T 



{1'J)J*+J + (1'J)*-1'(JJ*) 



1 e 4 



dr 



(2i'_j+ji + i'(jj*) ± - {i'j)±ji - 



Now we use that 



(73) 



l' + = -^£'(1 -cos0) = V2E'sm 2 ^-, l'_ = -^£'(l + cos#) = V2E'cos 2 °-. 
V2 2 y2 2 

We also direct qr along the x-axis, so that 

l' x = -E'sm9. 

The cross-section becomes 

da = -^—^V2E' cos 2 ^dr[(2J+j; + tan 2 ^(JJ*)± + v^tan °-{J x J* + + J+J*) 



(74) 



(75) 



1 e 9 9 9 

-^r—^V cos 2 -dr[(\V2J + + tan - J x \ 2 + tan 2 -|J„| 



The first factor in (76) is 



We present 



1 e 4 
2P +(? 4 



V2£' 



4^ 2 



cos 



2 " MV aMott - 



dH 1 

dT = J^f2E> dTh ' 



(76) 



(77) 



(78) 



where dr^ is the phase volume of the observed hadron. Then we finally get the exclusive 
cross-section as 



da = dE'dtt'dTh-^- (\V2J + + tan -J x \ 2 + tan 2 -I JJ 2 
AttM L vi 2 2 y . 



(79) 



To compare with the cross-section which is found by lifting the integrations in the inclusive 
cross-section we transform (79) to the system qx = 0. This is achieved by a rotation in xz 
plane by angle 9q determined by 



tan 6*o = — 



q x E' sin 9 
~q~ z ~ E - E< cos 9 ' 



(80) 
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We find 

where now 
with 



da = dE'dn'dT h 



&Mott 



AirM 



(|Z| 2 + tan 2 -|j/ 



1 

71 



1 ± cos T sin #o tan 



c + J + + c_ J_ + c x J; 




c x = sin #o + cos #o t an 



(81) 
(82) 
(83) 



On the other hand, starting from the inclusive cross-section (2), dropping the integration 
over drh and passing from the system q + = to the system = we obtain 



da = dE'dn'dT h 



&Mott 



(\z\ 



with 
where 

and 



AttM 

C+J+ + C-J- + C X X 



2tan 2 -\J y \ 2 



c± = -j={l ± cos^o), c x 



sm (po 



tan 0o = 



xM 

~Q~' 



One can trivially prove that 



cos <fio = cos #o — sin 6q tan -, sin 2 0o + tan 2 - 



sin 6q + cos #o tan - 



So it turns out that 



c± = c±, 



(84) 

(85) 
(86) 

(87) 

(88) 
(89) 



Comparing the correct cross-section (81) with the one (84) found from the inclusive cross- 
section, we see that they do not coincide, due to the difference in c x and c x . However if we 
integrate over the azimuthal directions of the observed nucleon, then both | J x \ 2 and | J y \ 2 go 
over to (l/2)|Jj_| 2 and the interference term with J+J* vanishes. Then both cross-sections 
(81) and (84) give the same result. 

Thus for the cross-sections averaged over the azimuthal directions of the observed nucleon 
one can use the inclusive cross-section with the dropped intergration over p\ z and p\±. How- 
ever if one wants to study also the azimuthal dependence, then one has to use the correct 
expression (81). 

9 Appendix 2. Some calculational details 

9.1 Singularity at 9 = 

We present the cross-section (55) in the form 



(2ir) 3 2p 10 da a\ 



d?pi 



d cos 9d(j> 
AAmE 2 J (l-cos#) 2 



X(9,cf>), 



(90) 
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where 

cos 2 {9/2) 



= \r o J ( 2^ ++ (g,pi) +tan 2 ^(g.pi) ) . (91) 



With a zero lepton mass, the integral (90) diverges logarithmically at 9 = (corresponding to 
Q 2 = 0, a real photon). With a non-zero lepton mass mi the divergence disappears, but the 
integrand still results strongly peaked at 9 = 0. To make its numerical calculation feasible 
one has naturally to separate the peak from the background. 
Presenting at small 6 

X(9, 4>) = ^ 2 ^i(0) = (1 - cos0)Xi(0) + (cos 9 - 1 + ^ X!(<f>), (92) 

we write the integral (90) in the form 

(27r) 3 2p 10 da _ a 2 em f r (l - cosflpfi^) X(9, 0) - (1 - cos 9)X l (<P) ] 

d* Pl ~4AmE 2 J acOSm n(l-cos9 + ^y + (1-cosfl) 2 \ ,( ) 

where 

_ mi ( 1 J_ 

In the first term the integrations over 9 and <f> can be done explicitly. The result cor- 
responds to the Weiszsaecker- Williams approximation. From our study it follows that the 
dependence of X\ (<f>) on the azimuthal angle is trivial 

X 1 ((f)) = cos 2 <f>X$ + sin 2 <f>X$ , (94) 

where terms with cos and sin come from W ++ and W yy respectively. So the integration over 
4> leads to the substitution 

X l (0) =- n(X$ + X$ ) = ttX 10 (95) 
Integrating over cos 9 we finally get 

(27r) 3 2p 10 Ar (2n) 3 2p 10 do- ww , a 2 em f X^j-Q^c^X^ 

<P Pl = d^l + 4A^J dCOs9d<P (1-cos^) 2 ' (96) 

where the Weizsaecker- Williams term is given by 

(27rf2 Pl0 da ww Tta 2 ev 



-X 1Q (\n^-iy (97) 



d 3 pi AAmE 2 \ /i 2 

The integration over 9 and 4> m the second term in (96) can only be done numerically. 

9.2 Singularities of the short-range part of X 

In the rescattering region the calculation of the short-range contribution to the integral (96) 
encounters difficulties generated by the singularities of Re at (q + p) 2 = m 2 . Denote 

X r (9, <j)) = X(9, <p)-(l- cos 9)X l (4>). 
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As a function of cj) X r has poles of the 1st and 2nd order in cos cj>: 

X r {0, </>) = X rr (9cf>) + — r f ] + ^ -2 , (98) 

A(9)-B(9)coscf> (a(0) - B(9)ccx<f>) 

where X rr is a regular function. The poles exist only in a certain interval of 9 

Gmin < $ < 6 max (99) 

determined by the condition 

\(f)\ 

< 1, (100) 



A(9) 



B(9) 

where rescattering is kinematically possible. 
Integration of Re X r over (f> evidently gives 

[ d(f>ReX r (9,(t>) = f dc/)ReX rr (9,(j)), (101) 

for 9 inside the rescattering interval (99), and 

/ # Re X,.(», *) = / #Re X A 0, 0) + 2* ( ^'.^ + ^ ( , )2 ! ( ^ )2 |3 /2 ) («*> 

for outside the rescattering interval. As we observe in the latter case the result is a non- 
integrable function of 9 at points 9 m i n and # ma;c . 

To overcome this difficulty we have to remember that in reality the double pole at A — 
B cos (f> is split into two close poles when one takes into account the small nuclear momenta 
neglected in in (102). Thus in the vicinity of, say, 9 m i n the integrand in (96) 

Y ^ = n / d ^ Re ' ( 103 ) 

(1 — cos 9y J 

will have a singularity structure 

Y{9) = Y r {9) + 1 ( - C - - C + Cl " ) , (104) 

1 \ V Qmin — 8 V #min + t(X — 9 J 

with Y r a regular function and a « 1 proportional to the difference (q, ki — ) . In the 
limit a^Owe obtain from (104) 

V Grain ~ & {Omin ~ 9) 3 / 2 ' 

where 

rj = a-c 1 , ( = ct/2, 

with a non-integrable singularity. However, correctly integrated, (104) has a finite limit at 
a — > due to the dependence on a of the limits of integration. Indeed one trivially finds in 
this limit 

/ Y{9)= Y r {9) + 2 v V9 min -9 - -=L=. (106) 
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Quite a similar formula is obtained for the integral over the interval Qraax 

< 9 < 6>i with the 

substitution 6 m in — 9q ^ 9\ — 9 max and of course different values for 77 and £• 

In principle the regularization functions ri 2(0) and numbers r/ and £ for the two singular 
points in 9 can be obtained analytically. However the resulting expressions are very cum- 
bersome, so that in practice we obtained them by studying numerically the corresponding 
singular behaviour. 

An additional difficulty consists in that the lower singular point 9 = 9 m { n for values of X\ 
close to unity turns out to be quite small and so very close to the the point 9 = at which 
function Y r aquires a singularity due to subtraction. To solve this problem we divide the 
interval [0, 9 m i n ] into two [0, 9$] + [9q, 9 m i n ] and apply the regularization (106) only in the 2nd 
interval. As a result the integrations over these two intervals give large numbers cancelling 
each other, which imposes very stringent requirements on the numerical precision. 
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11 Figure captions 

Fig. 1. Diagrams illustrating different mechanisms of cumulative production: spetator (a), 
direct (b), compressed tube plus rescattering (c and d). 

Fig. 2. Cumulative (solid line) and rescattering (dashed line) kinematical regions in the 
backward hemisphere at 6 GeV. 

Fig. 3. Effective internucleon momentum space potential as extracted from the photo- 
production data [ 9 ]. 

Fig. 4 Nuclear coefficients fo and cq. 

Fig. 5. Inclusive nucleon leptoproduction cross-sections on Cu at the incident energy 6 
GeV in the backward hemisphere. 

Fig. 6. Same as Fig. 5 on the Pb target. 

Fig. 7. Double differential cross-sections in pi and Q 2 for nucleon leptoproduction on Pb 
at the incident energy 6 GeV, emission angle 9\ = 90° and Q 2 = 1 (GeV/c) 2 . 
Fig. 8. Same as Fig. 8 at 0i = 180° and Q 2 = 1 (GeV/c) 2 . 
Fig. 9. Same as Fig. 8 at 9 X = 90° and Q 2 = 4 (GeV/c) 2 . 
Fig. 10. Same as Fig. 8 at 9 X = 180° and Q 2 = 4 (GeV/c) 2 . 
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Fig. 11. Double inclusive cross-section for the reaction e + 16 O — > e'(p e ) + N(p\) + X 
at E = 2.4 GeV, Q 2 = 0.8 (GeV/c) 2 and q = 0.439 GeV for the nucleons ejected into the 
backward hemisphere. Positive p m iss = |q — Pi | correspond to 4>\ = 180°, negative p m iss to 
(f>i = 0. The lower curve was obtained with an energy independent internucleon potential, 
the upper curve with a vector-like potential. Experimental points are from [ 11 ]. Lower 
(upper) points correspond to the lpi/2 (IP3/2) nuclear level. 
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